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DISCOVERY OF A GALAXY CLUSTER VIA WEAK LENSING 
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ABSTRACT 

We report the discovery of a cluster of galaxies via its weak gravitational lensing effect on background 
galaxies, the first spectroscopically confirmed cluster to be discovered through its gravitational effects 
rather than by its electromagnetic radiation. This fundamentally different selection mechanism promises 
to yield mass-selected, rather than baryon or photon-selected, samples of these important cosmological 
probes. We have confirmed this cluster with spectroscopic redshifts of fifteen members at z=0.276, with 
a velocity dispersion of 615 km s~^. We use the tangential shear as a function of source photometric 
redshift to estimate the lens redshift independently and find zi — 0.30 ± 0.08. The good agreement with 
the spectroscopy indicates that the redshift evolution of the mass function may be measurable from the 
imaging data alone in shear-selected surveys. 

Subject headings: gravitational lensing — galaxies: clusters: general 



1. INTRODUCTION 

Clusters of galaxies are essential tools for developing our 
understanding of structure formation and for probing cos- 
mological parameters. In particular, the redshift evolution 
of the cluster mass function is a sensitive diagnostic of f2„i , 
sufficiently sensitive that the existence of even one or two 
massive clusters at z ^ 0.8 favors a low-density universe 
(Donahue et al. 1998, Ebehng et al. 2000). This argu- 
ment assumes Gaussianity in the primordial fiuctuations; 
clusters are equally useful at constraining primordial non- 
Gaussianity given an independent measure of (Robin- 
son, Gawiser & Silk 1999). Precise estimates of either 
quantity will require large, unbiased samples of clusters at 
a range of redshifts. Locally, unbiased samples are cru- 
cial for measuring ag and for estimates of f2m that rely 
on the fair sample hypothesis, i.e. that the composition 
of clusters in terms of baryon fraction or mass-to-light ra- 
tio is representative of the rest of the universe (Evrard 
1997, Carlberg et al. 1996). Finally, the evolution of the 
numbers of cluster-sized masses as a function of redshift 
constrains the cosmological constant and the dark energy 
equation of state (Tegmark 2001). 

If clusters are to serve as measures of the clumping of 
mass, we must identify them directly from observations of 
the mass distribution. This is difficult, however, because 
the vast majority of their matter is dark. The traditional 
methods of discovering clusters rely on optical emission 
from galaxies [e.g. Abell, Corwin & Olowin 1989; Zarit- 
sky, et al. 1997) or X-ray emission from a hot intracluster 
plasma {e.g. Boehringer et al. 2000). Clusters with a lower 
mass-to-light ratio or more baryons may well be overrep- 
resented in these samples. Radiation from baryons is a 
complicated proxy for mass, and hydrodynamic baryon- 
CDM models have been proposed to study bias in clusters 
(Blanton et al. 1999). 



The development of weak gravitational lensing tech- 
niques (see Mellier 1999 and Bartelmann & Schneider 2001 
for reviews) suggests a different approach: shear selection 
(Schneider 1996). Because all types of matter participate 
in lensing, deep, wide-field imaging will reveal in its shear 
pattern all mass concentrations regardless of mass-to-light 
ratio or baryon fraction. Of course, no single technique 
will be completely unbiased. For example, shear selec- 
tion, like optical selection and the Sunyaev-Zel'dovich ef- 
fect (Joy & Carlstrom 2001), will tend to be more sensitive 
to line-of-sight superpositions of unrelated structures, to 
which X-ray selection, with its dependence on the square 
of the density, is less vulnerable. With X-ray and (until 
now) shear selection, optical foUowup is still necessary to 
provide a redshift. Comparison of differently selected sam- 
ples will always be necessary. Still, the baryon and photon 
independence provided by shear selection are powerful fea- 
tures which may produce surprising new samples. 

Two "dark" mass concentrations, found through weak 
lensing analyses of apparently unrelated previously known 
clusters, have been reported (Erben et al. 2000; Umetsu & 
Futamase 2000). Due to the absence of associated galaxies, 
the redshift and therefore masses of those clumps remain 
unknown. We report here the first spectroscopically con- 
firmed shear-selected cluster. We also introduce photomet- 
ric redshift techniques into the selection of sources for the 
weak lensing analysis. The recent development of these 
techniques (see Connolly et al. 1995; Hogg et al. 1998) 
greatly improves the promise of shear selection, as source 
galaxies may be divided into redshift bins to tune the sen- 
sitivity to lensing mass concentrations through a series 
of lens redshifts (mass tomography). In this paper, we 
use photometric redshifts to demonstrate that the shear- 
selected mass concentration is at roughly the same redshift 
as the cluster galaxies. This has always been an assump- 
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tion in lensing analyses of clusters, but this measurement 
implies that mass tomography is feasible. 

2. IMAGING OBSERVATIONS AND LENSING DETECTION 

We took BjVRI images of a "blank" field (containing 
no known cluster) centered at 23:47:46 +00:57:42 (J2000) 
using the Big Throughput Camera (BTC, Wittman et 
al. 1998) in photometric conditions at the Blanco 4-m 
telescope of Cerro Tololo Inter-American Observatory in 
1997 and 1998. Details of the data reduction, galaxy cat- 
alogs and seeing correction may be found in Wittman et 
al. (2000, WOO); here we give only the basic parameters. 
The final images are roughly 40' square with 0.43" sam- 
pling and 1.3" FWHM after point-spread function correc- 
tions, with object counts peaking at Bj = 26.4, V = 26.1, 
R = 25.6, and / = 24.4 (isophotal magnitudes). In the 
final catalog, each object's shape is a weighted mean of 
the shapes measured in the different filters, as described 
in WOO. 

In addition to the WOO magnitude cut of 23 < i? < 26, 
we imposed a color cut oi Bj — R < 1.5 to emphasize blue 
field galaxies at higher redshift, a tactic often used in weak 
lensing analyses of known clusters (Tyson, Wenk & Valdes 
1990). Using the method of Fischer & Tyson (1997), we 
constructed a mass map from the remaining 31,000 galax- 
ies (Figure 1). A mass concentration stands out near the 
southwest corner; its peak is significant at the 4.5cr level, 
based on mass maps of bootstrap resampled catalogs. The 
mass concentration disappeared, as it should, under two 
null tests: rotating each galaxy in the catalog by 45°, and 
assigning the shape of a random different galaxy to each 
galaxy position. To check for any bias casued by proximity 
to the field edge, we simulated a field of the same size and 
plate scale filled with random galaxies, distorted the image 
as if lensed by a cluster of moderate (800 km s~^) velocity 
dispersion, convolved and added noise to match the seeing 
and noise of the observations, cataloged, selected sources, 
and made mass maps as for the observations. Repeating 
this for a series of ten center-to-corner cluster positions, 
we found that the bias was small (< 10% in peak density) 
and in the sense of reducing, not enhancing, our sensitivity 
to mass concentrations near the corner. 

We then made a color composite image from the Bj 
and R images (Figure 1.5, separate jpeg). A concentra- 
tion of reddish elliptical galaxies appears near the posi- 
tion of the density peak (23:46:23.85 +00:45:00.8 for the 
brightest cluster galaxy versus 23:46:24.0 +00:43:58.8 for 
the density peak, a displacement of just over 1'). Nothing 
appears at that position in the ROSAT database, and the 
NASA/IPAC Extragalactic Database contains only one 
object in the area, a large spiral which appears to be a 
foreground field galaxy. We concluded that this candidate 
cluster was worthy of spectroscopic followup. 

3. SPECTROSCOPIC FOLLOWUP 

We designed a single slitmask containing 26 objects be- 
lieved to be cluster members, selected by avoiding blue 
objects and choosing suitably bright yellow/red galaxies 
in the color composite image in the area around the puta- 
tive cluster. This slitmask was used with the Low Resolu- 
tion Imaging Spectrograph (Oke et al. 1995) on the Keck 
Telescope on the night of Nov. 23, 2000. A single 1200 



sec exposure with the 300 1/mm grating was obtained; the 
spectral resolution was 10 A with a 1 arc-sec wide slit, and 
the region from 4000 to 8000 A was covered for each object. 
Two of the objects proved to be Galactic M dwarfs, while 
17 are members or near outliers of a cluster at z=0.276. 
One is a foreground galaxy, and six are background galax- 
ies. We assume an instrumental contribution to the ve- 
locity dispersion of 100 km/sec in the rest frame. While 
there are two outliers among the 17 possible cluster mem- 
bers, both sigma-clipping and the biweight estimator of 
Beers, Flynn & Gebhardt (1990) yield a velocity disper- 
sion estimate (in the rest frame) of — 615 ± 150 km 



4. LENS REDSHIFT AND MASS ESTIMATES USING 
PHOTOMETRIC REDSHIFTS 

With the cluster redshift in hand, we need only the red- 
shift distribution of the lensed sources to derive the clus- 
ter mass from the shear. In the past it has been diffi- 
cult enough to estimate the mean of this distribution, but 
photometric redshifts can provide a detailed distribution, 
in principle even an appropriate weight for each source 
galaxy. In this paper, we use tangential shear as a func- 
tion of source photometric redshift to estimate the lens 
redshift zi in a way independent of the spectroscopy. We 
demonstrate that the mass causing the shear signal is at 
roughly the same redshift as the cluster. 

For each galaxy observed in all four filters, we used 
the HyperZ package (Bolzonella, Miralles & Pello 2000) 
to compute a redshift probability density function (PDF) . 
We multiplied this by another PDF computed from the 
galaxy's apparent magnitude, assuming that a Schcchter 
(1976) luminosity function with M^. = —19.73 and a = 
— 1.28 (Folkes et al. 1999) holds at each redshift (we as- 
sume Hq = 70 km s-i Mpc"\ = 0.3, and f^A = 0.7 
throughout). The latter PDF is quite broad, but serves a 
valuable purpose by suppressing high-redshift peaks which 
often appear in the former PDF. The first and second mo- 
ments of the final PDF product provide the estimated red- 
shift and its statistical error. Henceforth, we use a catalog 
of 26,000 galaxies which is defined not by the magnitude 
and color cuts that went into Figure 1, but by the re- 
quirement of detection in each of four filters leading to 
a photometric redshift. Of these 26,000 galaxies, 13,000 
were also in the initial catalog. The median Zphot in this 
catalog is 0.58. 

We verified the accuracy of the photometric redshifts by 
comparison with spectroscopic redshifts of 31 galaxies in 
the range 0.23 < ^spec < 0.83. This sample was conrposed 
of 19 galaxies with redshifts described in Section || (the 
ones detected in all four filters), plus another 12 galaxies 
(0.24 < z < 0.83) in a different region of the 40' field, 
kindly provided by R. Guhathakurta. We find little bias, 
with (zspoc-^phot)/ (l+^^spcc) = -0.027 on average, and an 
rms of only 0.059 in the same quantity. A detailed analysis 
of this photometric redshift method and its application to 
larger datasets is in preparation (Margoniner et al. 2001). 

Before examining the dependence of the shear on source 
redshift, we must first account for the equally large de- 
pendence on projected position relative to the lens. We 
are not interested in a detailed reconstruction of the lens; 
rather, for a given Zphot we would like to collapse in- 
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formation from all sources at a wide range of projected 
radii into a single number characterizing the lensing sig- 
nal at that source redshift. Hence at a given Zphot we fit 
a very simple model, a singular isothermal sphere (SIS), 
to the radial dependence of the tangential shear (cen- 
tered on the peak in the mass map), and use the fit- 
ted amplitude and its uncertainty. The assumption of 
a particular profile should not introduce any bias as a 
function of 2;phot- To test this, we also considered NFW 
(Navarro, Frenk & White 1997) profiles. Although our 
data cannot constrain the scale radius Vg , for fixed (225 
kpc, Bartelmann, King & Schneider 2001), the results do 
not change significantly from the SIS case. We therefore 
choose the SIS for simplicity. For an SIS, the tangential 
shear 7t(r) = S(r)/Ecrit, where E(r) ^ r^^ is the pro- 
jected mass density, Scrit = {c^ / ^t^G){Ds/ DigDi) is the 
critical density, and Dg, Dig, and Di are the angular di- 
ameter distances from observer to source, lens to source, 
and observer to lens, respectively. Since only Ds and Dig 
are changing with source redshift, the amplitude of an 
fit to the tangential shear should grow with the distance 
ratio Dg/Dis. 

Figure 2 shows the fitted tangential shear (at a fiducial 
radius of 1 Mpc) as a function of Zphot- It is consistent 
with zero for Zphot < 0.3 and increases monotonically for 
0.3 < 2phot < 1-3 (the upper limit is to avoid extrapolat- 
ing too far beyond the range of spectroscopic verification 
of ^phoi)- The dotted line illustrates the shear expected 
from a lens at zi — 0.276, with S fit to the points. This 
is a good fit (xt ~ 1)- We explore the range of zi al- 
lowed by the shear data by stepping zi through the range 
0.025 < z; < 1.3 in steps of 0.025, and refitting at each 
step. The probabilities corresponding to s-t each step 
are plotted in Figure 3. The median and mode of this 
distribution are at zi = 0.31 and 0.30 respectively, with a 
68% confidence interval 0.225 < zi < 0.375 (these numbers 
change by less than 0.01 when an NFW profile is used). 
The fit for z; = 0.3 is also shown in Figure 2 (dashed 
line). Thus the lens roughly coincides with the cluster of 
galaxies in redshift as well as in right ascension and decli- 
nation. This method can be used to estimate the redshift 
of any newly discovered lensing mass from the lensing data 
alone. Such a procedure may well become a standard part 
of shear-selected cluster surveys. 

Finally, we estimate the mass and mass-to- light (M/L) 
ratio using the best-fit projected mass. Still assum- 
ing an isothermal sphere, the projected mass inside of 
250 kpc (where it is convenient to measure the light) is 
2.8 ± 0.6 X IO^'^Mq, assuming zi = 0.276. For the range 
of zi allowed by the 7t(2;phot) curve, Mpi.oj(< 250 kpc) 
ranges from 1.8 — 3.7 x lO^^M©. The velocity dispersion 
implies Mproj(< 250 kpc) = 0.7 ±0.3 x lO^-^M© under the 
SIS assumption. The discrepancy may be due to the SIS 
assumption: Unlike the redshift estimate, the mass esti- 
mate is sensitive to the profile assumed. Converting from 
observed / band to rest-frame R using the approach of 
Fischer & Tyson (1997), we find M/Lr = (560 ± 200)h 



in this region. Compared to other clusters (Mellier 1999), 
this is high but not exceptional. 

5. DISCUSSION 

Baryon-unbiascd samples of mass concentrations over 
a wide range of redshift will be of critical importance in 
constraining cosmological parameters. To realize the po- 
tential of this technique, weak lensing observations must 
have the sensitivity to discover clusters over a broad part 
of the cluster mass function. For example, constraints on 
^Im from the mass function at high redshift currently in- 
volve only a few very massive clusters, and such extreme 
clusters lie far out on the tail of the mass distribution, 
which may not be Gaussian. Weak lensing surveys can 
probe clusters an order of magnitude less massive. Un- 
like X-ray and optical selection, a shear signal does not 
diminish as the square of the luminosity distance, so that 
low-mass clusters should be detectable even at high red- 
shift as long as the photometric redshifts are accurate in 
eliminating foreground sources. 

We have demonstrated the serendipitous discovery, with 
high signal-to-noise, of a rather modest cluster via a weak 
lensing analysis of a single 40' field. The cluster was spec- 
troscopically confirmed at z = 0.276, with a velocity dis- 
persion of 615 km sec~^. The tangential shear follows 
source (photometric) redshift in a manner which requires 
the lens to lie at or near the cluster redshift. Thus, all 
the ingredients are in place for a truly shear-selected sam- 
ple of clusters, in which any putative mass clump can be 
confirmed, and its redshift estimated, from the multi-color 
imaging data alone. This technique is also capable of an- 
swering the question of the existence of any truly dark 
clumps. The redshift (and therefore mass, and mass-to- 
light ratio lower limit) of any such clumps can only be 
derived from the shear versus source redshift curve. 

This further suggests that the promise of three- 
dimensional mass tomography (Tyson 1995, 2000) over 
cosmologically significant volumes can be realized in wide- 
field deep imaging surveys. Note that such a cluster is 
not unexpected in the volume probed by this data (Rah- 
man & Shandarin 2001). Ongoing cosmic shear surveys 
covering tens of square degrees {e.g. the Deep Lens Sur- 
vey^, DESCART^) should discover significant samples of 
shear-selected clusters (Kruse & Schneider 1999) and be- 
gin to constrain flm and dark energy through the redshift 
evolution of the cluster mass function. 

We thank P. Guhathakurta for providing spectroscopic 
redshifts outside the cluster, and the CTIO staff for their 
help with the BTC project and for their upgrading and 
maintenance of the delivered image quality of the Blanco 
telescope. Cerro Tololo Inter- American Observatory is 
a division of National Optical Astronomy Observatory 
(NOAO), which is operated by the Association of Univer- 
sities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. BTC 
construction was partially funded by the NSF. 
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Fig. 1. — Projected mass map, smoothed on a 2' scale, of the 40' field. Higher-density regions are shown darker. Contours 
are equally spaced in arbitary units (but linear in projected mass density); negative contours are omitted for clarity. One 
peak, at lower right, stands out at twice the density of any other peak (4.5fT). This mass overdensity corresponds to a 
small cluster of galaxies spectroscopically confirmed at z=0.276. The width of this field at that redshift is 13 Mpc. North 
is up and east is to the left. 
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Fig. 2. — Measured shear as a function of source photometric redshift (points). Each point represents the amphtude of 

a best-fit isothermal shear profile at a fiducial radius of 1 Mpc, with vertical error bars indicating the uncertainty in the 
fit. The dotted and dashed lines represent the shear expected from lenses at z = 0.276 and z — 0.30, the spectroscopic 
and shear-derived redshifts respectively (the different amplitudes reflect slightly different best-fit masses). The horizontal 
error bars represent the nominal widths of the bins only; the effect of scatter in the photometric redshifts is neglected. 
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Fig. 3. — Lons rcdshift probability density function (unnormalizcd) . Solid line; at each putative lens rcdshift zi, we fit a 
lens mass to 7t(2:phot) (Figure 2) and compute the probability from the for the remaining four degrees of freedom (five 
data points minus one fit parameter). Dotted line: we repeated the process assuming an NFW profile with = 225 kpc. 
Either assumption leads to a most probable zi within 0.03 of the spectroscopic value. 
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